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Abstract 

The purpose of this paper is to describe asymptotic formulas for determinants 
of certain operators that are analogues of Wiener-Hopf operators. The determinant 
formulas yield information about the distribution functions for certain random variables 
that arise in random matrix theory when one rescales at "the edge of the spectrum" . 

1 Introduction 

This paper is concerned with the asymptotics of Fredholm determinants of operators that 
arise naturally in random matrix theory and are similar in many ways to finite Wiener-Hopf 
operators. The operators, denoted by A a (f), are integral operators on £ 2 (R) with kernel 
given by 

f(x/a) / A(x + z)A(z + y)dz (1) 



where 

A{x) = / e it /3 e itx dt, 

and / G L°°(R). The function A(x) is the Airy function, generally denoted by Ai(x), and 
for this reason we call our operators A a (f) Airy operators. We will refer to the function / 
as the symbol of the Airy operator. 

If the term / °° A(x + z)A(y + z)dz in (0) is replaced by the sine kernel 

sin7r(x — y) 
x-y 
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the resulting operator has the same Fredholm determinant as a finite Wiener-Hopf operator, 
whose asymptotics are very well known ||. The similarities with the Wiener-Hopf (i.e., sine 
kernel) case become less surprising after the observation that 

f°° *, \ Al \ . A(x)A'(y) - A(y)A'(x) 

/ A(x + z)A(y + z)dz = — — — — — . (2) 

Jo x — y 

The proof of this as well as many other facts about the above kernel can be found in |5|. 
We shall use the fact that A(x) is rapidly decreasing at +00 and 0(|x| -1 / 4 ) at —00. For the 
complete asymptotics of the Airy function, we refer the reader to H. 

Under appropriate conditions A a (f) is a trace class operator, and thus the Fredholm 
determinant 

det(/ + A a (/)) 

is defined. The main goal of the paper is to compute the asymptotics of this determinant as 
a — > 00. 

The motivation for finding an asymptotic formula for the Fredholm determinant comes 
from random matrix theory, in studying so-called linear statistics, which are certain functions 
of the eigenvalues of random matrices. After a rescaling at "the edge of the spectrum" their 
characteristic functions become Fredholm determinants of our Airy operators. For general 
information about random matrices, we refer the reader to 0. For information about the 
connection of random matrices, characteristic functions and the Airy operators we refer the 
reader to || and ||. 

The paper is organized as follows. In the second section we derive the basic properties of 
A a (f) and related operators. In the third section we prove through a series of lemmas that 
for appropriate functions / and F 

Jim tr [F(A a (f)) - A a (F o /)] = tr [F(W(g)) - W(F o g)\, 

where W(g) is the Wiener-Hopf operator with symbol g(x) = f(—x 2 ). (The precise definition 
of W(g) will be given at the end of the next section.) The trace of the second operator on 
the left is easy to compute asymptotically. Taking F(z) = log(l + z) and using the known 
formula for the trace on the right side, we find that the Fredholm determinant is given 
asymptotically as a — > 00 by 

det (J + A a {f)) = exp {ci a i/2 + c 2 + o(l)} , (3) 

where 

1 f°° 

ci — — \fx log(l + f(—x)) dx, 

7T JO 

C2 = - / x (G(x)) 2 dx, 
2 Jo 
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and 



1 r 

G <*> = 2?/- 



■oo 



log(l + f(-y 2 ))dy. 



— oo 



This is proved under the assumption that / is a Schwartz function (although we could get 
by with much less) and 1 + f(x) 7^ for x < 0. 

This formula bears a strong resemblance to the corresponding asymptotic formula in the 
classical finite Wiener-Hopf case. The most notable difference is the power o?l 2 in the first 
term of the asymptotics. 

In the last section we describe the implications of formula @ for random matrix theory. 
The formula, as in the analogous Wiener-Hopf or Bessel kernel case (see H for details), 
proves that the distribution functions for certain linear statistics, now scaled at the edge 
of the spectrum, are asymptotically Gaussian. The recurrence of the Gaussian distribution 
highlights the universality seen again and again in random matrix models. 

2 Basic properties of the Airy operator 

We begin by defining the Airy transform A. For g 6 L 2 (TVj we define A(g) by the formula 



T~ x is the inverse transform and denotes multiplication by the function h(t) = e** 3//3 . We 
will also use the standard notation g and g for the Fourier transform and inverse transforms 
respectively. Observe that for g G L 1 n L 2 we have A(g)(x) = /f^ A(x + y)g(y)dy and, just 
as in the Fourier transform case, A(g) is the L 2 limit of J^ B A(x + y)g(y)dy as B — > 00 for 
all g e L 2 . 

Lemma 2.1 The Airy transform is unitary on L 2 and satisfies A^ 1 = A. 



Proof. Clearly A is unitary since T and Mh are. The Fourier inversion formula says that 
T~ x = J T ' = T 3 where Jg{x) = g(—x). It follows that 



Given the above definition of the Airy transform, we see that the Airy operator defined by 
(0) is alternatively defined as the operator Mf a APA, where P is multiplication by Xr+ and 
fa{x) — f{ x / a )- K is this representation of the operator which we will use. For appropriate 
/ this operator in turn will have the same Fredholm determinant as AMf a AP. We next 
derive a representation for the kernel of AMf a A for a class of functions /. 





A = T J M h J T = TM h -x T = A~\ 



□ 
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Lemma 2.2 Suppose that f is the inverse Fourier transform of a finite measure fi, 

f( x ) = — e^dfx(C). 



2tT J-oo 

Then the kernel of the operator AMfA is given by the formula 

-i(x+y)£/2 J(x-y) 2 /^ 



y/8lT J-oo y/W+O 



e -*x+uwi e ^-yn^dfj,^), (4) 



where ^Ji^ + is defined by taking arg(z£ + 0) equal to tc/2 when £ > and —n/2 when 
£<0. 

Proof. Consider first the case where /i is a unit mass at the point rj, and let K v be the 
operator with the corresponding kernel (|4]). Then f(x) is the function e v (x) = e lvx /y2ir. By 
Lemma 2.1 we see that we have to show 

K v = T~ x M h T~ x M e „ T M h -i T, 

or equivalently 

Let us compute both sides applied to a function in L 2 . 

Notice first that JF _1 M e T takes a function into + r))/y/2n. Therefore 

Mh T~ x M £rj T M h -i takes <p(£) into 

1 e ^ 3 /3 e -,(^) 3 /3 ^ + ^ = _L e -*(^+^ 2 +W3) ^ + ^ 



V27T V27T 

Hence 

1 /"OO 



T~ x M h T~ x M eri FM h -, ip{x) = — f e** e ~^+^W/V ^ + ^ d {= 

27T J-oo 
-i(r?x+r? 3 /3) .00 



27T 

On the other hand, we have 



e**e-* 6i r-W(p{£)dZ. (5) 



i „ — ir? 3 /12 i /.oo 

n V J v^vrv^TO v^f 7-00 ^ 7-oo y 

(If ip is a Schwartz function, for example, the interchange of order of integration this involves 
can be justified by integration by parts, and it suffices to show our two operators agree when 
applied to Schwartz functions.) The inner integral is easily computed and found to equal 



Thus we see that K^T^ 1 is equal to the right side of ©. 

The lemma is established for the special case of a unit point mass, and so for any linear 
combination of these. To establish the general result we approximate our given /i by a se- 
quence {/x n } of linear combinations of point masses such that the measures /i n are uniformly 
bounded and / </?(£) cf/i n (£) — > f </?(£) d[i(£) for any function ip which is bounded and con- 
tinuous. Then the corresponding functions /„ converge boundedly and pointwise to / and 
so the operators AMf n A converge strongly to AMfA. For the corresponding operators K n , 
it is easy to see that for Schwartz functions g± and g 2 we have (K n gx, g 2 ) — > (K gi, g%), so 
K n —>■ K weakly. Hence, since AMf n A = K n for each n, we have AMfA = K. □ 

To end this section we recall the definition of a Wiener-Hopf operator and certain of its 
properties. For a function g e L°°(R) the operator W(g) on L 2 (R + ) (which we identify with 
the functions in L 2 (R) which vanish on R~) is defined by 

W(g) = PT~ x M a TP. 

This is the Wiener-Hopf operator with symbol g. Notice the analogy with the operators 
PAMfAP. The fact that there is more than just an analogy will become apparent in the 
next section. One often sees a Wiener-Hopf operator defined as an operator on L 2 (R + ) 
with kernel of the form k(x — y) where k G L 1 (R). This operator is equal to W(g) with 
g(x) = J e~ lxu k{u) du. 

We state as a lemma two basic facts about Wiener-Hopf operators. 

Lemma 2.3 a) The spectrum of W(g) is contained in the convex hull of the essential 
range of g. 

b) If g is continuous and g(±oo) = then A is not in the spectrum ofW(g) if and only 
if \ ^ 0, g — A 7^ and 

^9-^)-=^- A arg (g(x) - A) = 0. 

Z7T — oo<x<oo 

3 Trace norm estimates and the Airy limit theorem 

We assume from now on that / is a Schwartz function. The reader can verify that this 
requirement is too restrictive and can, for example, be replaced by a weighted space condition. 
However, assuming that / is a Schwartz function simplifies the proofs and increases the clarity 
of the arguments. 

Recall that the Airy operator A a (f) equals Mf a APA and is thus similar to AMf a AP, 
which in turn is to unitarily equivalent to the operator U~ l AM f a APU where U is the unitary 
operator defined by Ug(x) = a l l 4 g(x^fa). Note that U commutes with P. This operator will 
act as a replacement for the Airy operator in the final computations. The next lemma, which 
involves a modification of the above operator, will be important in those computations. 
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Lemma 3.1 The operator U 1 (I—P)AMf a APU converges in the trace norm to the operator 
with kernel 



V8n J-oo + 

as a — > oo. 



Proof. By Lemma 2.2 the kernel of the operator (7 — P)AMf a AP is given by the formula 

J-oo \/2c + U 



and thus the kernel of £/ -1 (7 - P)AM fa APU is given by 

V8n J-oo v^TO 
Changing x to —x in the kernel for convenience gives 



j-oo \A£ + 



1 A 00 .^ILe-^/^e^-^/^e^^dexR+^XR+Cy). 

V^TT i-oo V^TO 

We shall show that replacing by 1 each of the two factors in the integrand which involve 
a leads to an error which is the kernel of an operator having trace norm o(l). We will use 
the general fact that the trace norm of an operator with kernel K(x,y), where y is confined 
to a set J, is at most || K || 2 + \J\ || dK/dy || 2 , where the norms are Hilbert-Schmidt norms. 

We first look at the error kernel arising from the replacement e - l ( x -y)£/ 2a3/3 — > 1 ; which 

is 

-W- r -Jjfi=e-^l^ ( e ^/ 2 " 3/2 -l)e^)V^^ XR+ ( x ) XR+(y ). 

V87T J-oo + 

This we call K(x,y) and find bounds on K and dK/dy. 

Clearly K(x,y) = 0(\x — y\/a 3 ^ 2 ) = 0(w/a 3 ^ 2 ), where w = x + y. We use this estimate 
for w < 1. To get a better estimate for w > 1 we write the kernel as a constant times 

J_ Z" 00 /(Q p -^/12a 3 / i(x-j/)C/2a 3 / 2 _ ^ ,2 ^_ i(x+y) 2 /4£ .£ 

w 2 J-oo v^FFo v J ^ d£ s 

and integrate by parts to obtain a constant times 



Of course we apply the product rule. Differentiating the various factors in the brackets 
leads to an extra factor a -3 or wa~ 3 ^ 2 , aside from the external factor 1/w 2 . If w > 1 we 
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see therefore that integration by parts yields a factor a -3 / 2 < w >~ 1= a~ 3 / 2 (l + w 2 )^ 1 ^ 2 
assuming that of course / is in a Schwartz function. Integrating by parts once more leads 
to a factor of a~ 3 < w >~ 2 . Thus for all positive x and y (w > 1 or w < 1) we see that our 
kernel satisfies K(x,y) = 0(a~ 3 / 2 < w >~ 2 ). 

We also have to estimate dK(x,y)/dy. If we differentiate 



we are left with two integrals. In one integral we get an extra ^/a 3 ^ 2 in the integrand, 
and the factor e l( - x ~ y ^^ 2a3 — 1 is replaced by e l ^ x ~y)^/ 2aZ 2 . As before this can be seen to 
be 0(a~ 3 / 2 < w >~ 2 ). In the other integral we get an extra w/£ which changes the factor 

. - ( . 

Here after three integration by parts we arrive at an estimate of 0(a~ 3 / 2 < w >~ 2 ). 



We have shown that 



d ( a- 3 ' 2 

K(x,y), —K(x, y ) = 0' 



dy ' \l + x 2 + y 2 J 

If we use the general trace norm estimate stated above, taking J = (k, k + 1) for k — 0, 1, • • • 
and adding, we find that the trace norm of the error operator is 0(a~ 3 / 2 ). 

If we consider the error due to the replacement e~ %( - l V2a — > 1 the argument is essentially 
the same and we find a bound for the trace norm of the resulting kernel of 0(a~ 3 ). This 
completes the proof. □ 

Here and below we shall use the notations Oi(-) resp. o±(-) to denote families of operators 
depending on the parameter a whose trace norms are O(-) resp. o(-). 

Lemma 3.2 . We have PAM fa AP = Oi(a 3/2 ) in general and PAM fa AP = oi(l) if f 
vanishes on R~. 

Proof. The kernel of our operator on L 2 (R + ) equals 

f(z/a) A(x + z) A(y + z) dz. 



I 

j — i 



For fixed z the kernel f(z/a) A(x + z) A(y + z) is a separable rank one kernel. To compute 
its trace norm, observe that by the estimates on the Airy function we have 



/ A(x + z) 2 dx 
Jo 



' 0(e- z ), z>0, 
k 0{< z > 1 / 2 ), z < 0. 
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It follows that the trace norm of our operator is at most a constant times 



J —oo 



f \f(z/a)\ < z >^ 2 dz + 



/■oo 

/ \f(z/a)\e-dz, 
Jo 



and the assertions of the lemma follow easily. 



□ 



Corollary 3.3 A a (f) is a trace class operator. 

Proof. After the replacement x — > —x the kernel in the statement of Lemma 3.1 becomes a 
Hankel operator with smooth kernel and thus is well-known to be trace class. The lemma 
implies that — P)AMf a APU is trace class. Thus (/ — P)AMf a AP is trace class. 

Lemma 3.2 tells us that PAMf a AP is trace class. Hence so is AMf a AP, and A a (f) = 
Mf a APA is unitarily equivalent to this. □ 

We remark that this argument could have been made much earlier. However it would 
have involved the same sort of estimates as in the proof of Lemma 3.1 and there was no 
point in doing this twice. 

The operator with kernel 



is a convolution operator, its kernel is a function of x — y. We will denote it by Kf. Thus 
PKfP is a Wiener-Hopf operator with symbol 



Thus PKfP = W(g), and the connection to Wiener-Hopf operators is now apparent. 

Lemma 3.1 told us that the operator — P)AMf a APU converges in the trace norm 

to (/ — P)W(g)P. The next lemma concerns the strong convergence of the operator 





Now 




and therefore our symbol is given by 




B a (f) = U- x AM ja AU. 



This is the last technical lemma before we can put the pieces together. 
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Lemma 3.4 The operator B a (f) converges strongly to Kf as a — > oo. 
Proof. We have to show that for any ip e L 2 (R) 



oo roo 



oo •/ —oo 



< & ! e -^+y^ 2a3/2 <p(y) dydti 



r r -01= <p{y) dy d£. 

J-oo J-oo +0 



in L 2 (R). We can restrict ourselves to a dense subset of ips since the B a (f) have uniformly 
bounded norms. 

Write the double integral on the left as 

+ J-oo 

This minus its purported L 2 limit equals the sum of the two error integrals 



7-oo 



and 

/(£) 



, ( e -^/12« 3 e - 4 x ?/2a 3 /2 _ x) f°° e Ha,-y)*/4t / x dyd £ / 7) 

y/i^ + J-oo 

The operator with kernel e l ^ x ~ y ^ 2 ^ is unitarily equivalent to, and therefore has the same 
norm as, the operator with kernel |£| 1//2 e %<yX ~ y f . Thus it has norm 0(|£| 1/2 ). The function 

(e-^ 3/2 - 1) <p{y) 

has norm 0(£/a 3//2 ), assuming as we may that y (p(y) G L 2 , and it follows that the inner 
integral in (H) has norm 0(|£| 3//2 /a 3 / 2 ). Hence (|]) itself has norm at most 0(l/a 3 ^ 2 ). 

As for ([7J), the inner integral equals a function ip^(x) whose L? norm is Od^) 1 / 2 ). Write 
(0) as the sum 

/(£) /-i£ 3 /12a 3 T^-ii^a 3 / 2 ., , f°° /(£) 



The L 2 norm of the function ( e -^ 3 / 12a3 _ i) e -^/2« 3/2 ^(x) is 0(|£| 7 / 2 /a 3 ) and so the first 
integral in (|7]) is a function whose norm is 0(l/a 3 ). As for the second integral, observe that 

|| ( e -**/«_i)^( x ) || 



9 



is 0(£) for all £ and a and tends to as a — > oo for each £ (by the dominated convergence 
theorem). Hence the integral obtained by taking norm under the integral sign in the second 
integral tends to as a — > oo, again by the dominated convergence theorem. This establishes 
the claimed strong convergence. □ 

Now we are ready to begin the final steps in proving (|3|). The operator (A a (f)) n = 
(Mf a APA) n has the same trace as (PAMf a AP) n which in turn has the same trace as 
(PU^AM^AUP) 71 = (PB a (f)P) n . In fact for any analytic function F defined on the 
neighborhood of the spectra of both operators and satisfying F(0) = we have 

tr F(A a (f))=tvF(PB a (f)P). 

It is the asymptotics of this last trace we shall compute. We think of our operators as acting 
onL 2 (R+). 

In the following two lemmas A will be in the resolvent set of the Wiener-Hopf operator 
PKfP = W(g). By Lemma 2.3b this implies that A is not in the range of g(x) = f(—x 2 ), 
so / _1 ({A}) is a compact subset of (0, oo). We can find a Schwartz function / which never 
takes the value A and which equals / outside some larger compact subset of (0, oo), and 
we can find one / which serves for all A in any given compact subset of the resolvent set of 
W(g). This will be our notation in what follows. 

Lemma 3.5 Let A be in the resolvent set of PKfP = W(g). Then A is also in the resolvent 
set of PB a (f)P for sufficiently large a and the inverses have uniformly bounded norms for 
A lying in any given compact subset of the resolvent set. 

Proof. Observe that for Schwartz functions fx and f 2 , since B a (fi) B a (f 2 ) = B a (fi / 2 ), 
PB a {h)PB a {f 2 )BP - PB a {h f 2 )BP = PB a {h){P - I)B a {f 2 )BP 

-> PK h {P - I)K h P = PK h PK h P - PK hh P = W( 9l ) W(g 2 ) - W( 9l g 2 ) (9) 

in trace norm since (J — P)B a (f 2 )BP — > (I — P)Kf 2 P in trace norm by Lemma |3.1| and 
PB a (f\) — > PKfc strongly by Lemma |3.4[ This also holds if the fa are constants plus 
Schwartz functions. 

We take fx — f — A and f 2 = (/ — A) -1 . Observe that the u g" corresponding to f 2 is 
(g — A) -1 , so in this case the relation @ reads 

PB a (f - X)PB a ((f - \)- l )P - PB a ((f - A) (/ - A)- X )P - W(g - A) W((g - A)" 1 ) - P. 
Now 

PB Q ((f - A) (/ - A)- X )P = P + PB a ((f - A) (/ - A)- 1 -1)P = P + 0l (l), 
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by Lemma 3.2. We conclude that 

PB a (f - X)PB a ((f - A)- X )P = W(g - A) W((g - A)- 1 ) + 0l (l). (10) 

The analogous formula holds when / — A and (/ — A) -1 are interchanged. The Wiener-Hopf 
operators W(g — A) and W((g — A) -1 ) are invertible, the first by assumption and the second 
since ind (g — A) -1 = — ind (g — A) = 0. The norms of the inverses are bounded uniformly in 
a and A lying in a compact set and the operators PB a ((f — A) _1 )P are uniformly bounded. 
This completes the proof. □ 

Lemma 3.6 Suppose F is analytic in a neighborhood of the spectrum ofW(g). Then we 
have as a — > oo 

F(PB a (f)P) - PB a (F o f)p = F(W(g)) - W(F o g) + 0l (l). (11) 

Proof. By Lemma 3.5 PB a (f — X)P is invertible for sufficiently large a with uniformly 
bounded norm for A lying in a compact set in the resolvent set of W(g). Let A also be in the 
domain of F. Then 

(PB^-\)P)- 1 -PB a ((f-\)- 1 )P=(PB^-\)P)- 1 [I-W(g-\)W((g-\)- 1 )} + o 1 (l) 

= W(g - A)" 1 [I - W(g - A) W((g - A)" 1 )] + 0l (l). 

The first equality follows from ([10]) and the uniformity of the norms of the inverses. The 
second equality uses the strong convergence of (PB a (f — A)P) _1 to W(g — A) -1 and the fact 
that I - W(g - A) W{{g - A)" 1 ) is trace class. Thus 

(PB a (f - A)P)- X - PB a ((f - A)^)P = W(g - A)" 1 - W((g - A)" 1 ) + 0l (l). 

Multiplying by P(A) and integrating over an appropriate contour gives 

F(PB a (f)P) - PB a (F o f)P = F(W(g)) - W(F o g) + 0l (l) 

for any F analytic in a neighborhood of the spectrum of W{g). □ 



We will be interested in the trace of the first operator on the left side of (0). The next 
lemma will tell us the trace of the second operator. 

Lemma 3.7 For any Schwartz function f we have 

v 3/2 



tr PB a {f)P = / y/xf(-x) dx + o(l) 

7T JO 
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Proof. The kernel of PB a (f)P equals 

i r /(0 



V^TT J-oo 

and thus 



/■°° e -^ /12Q 3 e _ l(x+?;)?/2a 3/ 2 e ^_, m ^ Xr+ ( x ) Xr+ ( 2/ ) ; 

J-oo Jit + 



tr PB(f)P = — / / ffO e -^/i2o» e -^/ Q 3/ 2 ^ 

V87T JO J-oo + 



We write this as 

a 3/2 ^ ^ 



JO J-oo \ It + 



\f%1\ JO J-oo + 

and then replace the term e -*f 3 / 12 a 3 by ]_ j us ^ as j n Lemma |3.1| to find that the trace is given 
by 



We can write this in a more familiar form by replacing the term 1 / + in the above 
integral with Jf^ e _m ^du. Integrating over £ we find that this equals 

q/3/2 rCO roo 



1lX Jo J-oo 

or 



3/2 



f(—u 2 — x)du dx + o(l) 

/•oo 

/ v^r /(— x) dx + o(l) 

7T JO 

as claimed. □ 

We now derive our main result on determinants of Airy operators which gives the 
promised formula for the asymptotics. We assume, as always, that / is a Schwartz 
function. 

Theorem 3.8 Assume g(x) = f(—x 2 ) ^ —1. Then as a — > oo 

det (/ + A a (f)) = exp {ci « 3/2 + c 2 + o(l)} , (12) 

1 f°° 

ci = - / l°g(l + f(~ x )) dx, 

IX Jo 

I f°° / \ 2 

° 2 = 2Jo X \^° g(yl + 9))( x )j dx - 
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Proof. Assume first that \\g\\oo < 1- Then by Lemma 2.3a the spectrum of W(g) lies in 
the open unit disc with center 0. Therefore F(z) = log (1 + z) (the branch equal to when 
z = 0) is analytic on the spectrum and so we may apply Lemma 3.6. This and Lemma 3.7 
tell us that there is an asymptotic formula of the form (|L2|) where C\ is as stated and 

c 2 = tr [log (/ + W(g)) - W(\og (1 + <?))]. 



It is known that this equals the expression given for c 2 in the statement of the theorem [|10 

To remove the restriction on g we introduce a parameter t and we would like to define 
a family of functions ft by 1 + ft — e tl ° z<yl+ ^\ so that for small enough t our asymptotic 
formula holds. The problem is that —1 may lie in the range if /, and even if it didn't we 
might not have i(l + /) = 0, which is what we need to define a logarithm which is a Schwartz 
function. So, as in the preceding lemmas, we introduce a function / which equals / except 
on a compact subset of (0, <x>) such that 1 + / ^ and ind (1 + /) = 0. Then we define f t 
by 

1 + /t = e tiog(i+/) + / _ / ; 
Of course fi — f- Moreover, with g t (x) = ft(—x 2 ^ 



I ■ 



l+g t = e t log(1+9) 

for all t. For sufficiently small t we have H^Hoc < 1 so that our formulas hold. 

Observe that det (I + A a (f t )) is a family of entire functions of t depending on the pa- 
rameter a. Suppose we can show that 

det(I + A a (ft)) = 0(e tc ^ /2 ) (13) 
for large a uniformly on compact t-sets. Then the limit relation 

Jim e-' ClQ3/2 det (I + A Q (ft)) = e tc \ 

which we know holds for sufficiently small t, will hold for all t and therefore t = 1. 

To prove (O) we go back to the PB a (f)P regarded as operators on L 2 (R + ). We have 



det (/ + A a (f t )) = det PB a {\ + f t )P. 

Now i(l + g t ) — and so W(l + gt) is invertible by Lemma 2.3b. Therefore by Lemma 3.5 
with F(z) = z~ l we know that PB a (l + f t )P will be invertible if a is large enough. (This 
will hold for all t in any given compact set.) For these a we have 

j t log det (I+A a (ft)) = tr \{PB a {l+f t )Py l £pB a (l+f t )P} = tr \{PB a {l+f t )P)- 1 PB a (h t )P], 
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where 

ht = log (1 + /) e ' log{1+/) . 
By Lemma 3.6 with F(z) = z~ l we know that 

(PB a (l + f^P)- 1 = PB a ((l + f^P + 1 (1). 

Also, by (|), 

PB a ((l + f t y l ))PB a {h t )P = PB a ((l + f^- 1 ht)P + O x (l), 
so that we have shown 

j t log det (I + A a (f t )) = tr PB a ((l + ft)' 1 h t )P + 0(1). 

But Lemma 3.7 tells us that with an error o(l) 

tr PB a (l + f^ 1 ht)P = / {l + M-x^hti-x) ^dx 

IT JO 

/ \fx log (1 + f(—x)) dx. 
Jo 



TC JO 

because / = / on Rr. Thus, 

- log det (J + A a (ft)) = / log (1 + f(-x)) dx + 0(1). 

at 7T JO 

Integrating over t from to t and exponentiating gives (p~3|) and completes the proof. □ 

4 Applications to random matrices 

Theorem 3.7 can be applied to find limiting distribution functions for a class of random 
variables which are functions of the eigenvalues of a random matrix. In many different 
ensembles of matrices it has been shown that the distribution functions are asymptotically 
normal [|I|, |2|, |3], and this will be shown also to be the case in our examples. The term 
ensemble refers to the probability density assigned to some space of matrices, and this in 
turn induces a density on the space of eigenvalues of the matrices. For the Gaussian Unitary 
Ensemble (GUE) the density on the space of eigenvalues is given by 

Pn{xi, x N ) = -^y det K(xi, Xj) \^ j=1 (14) 
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where 

N-l 

K N (x,y) = <l>i( x )Mv)- ( 15 ) 

and <fti is obtained by orthonormalizing the sequence {x l e~ x2 ^ 2 } over R. If N is large it is also 
well known that the density of the eigenvalues is supported on approximately the interval 
(-V2N, V2N). These facts can be found in 0. 

The random variables of interest here are ones that are often called linear statistics and 
are of the form 

N 

E/(V«). 
i=i 

where A, are the eigenvalues and / is an appropriate function. Our goal is to study these 
random variables applied to the eigenvalues near the edge of the spectrum and to this end 
we rescale and replace the sum by 

N 

J2f(2 1/2 N^(\ t -V2N)/a). 
i=i 

The purpose of the translation by the term y2N is to move to the edge and the factor 
2 1 / 2 jV 1 / 6 has the effect of making the eigenvalue density of the order 1. Otherwise the 
eigenvalues "bunch up" or "spread out" and all the results become more or less trivial. 

To study the distribution function of this random variable we use its characteristic func- 
tion, or inverse Fourier transform. This characteristic function is given by 



>*(*)=/ •■•/ e-^-^^-^^P^,...^^^---^ 



It is a general fact that 

/OO POO N 

" " / IK 1 + 9(xi)) P N {xi, ...,x N )dx x --- dx N = det(I + gK 
-OO J — CO ■ i 



N) 

i=l 

where g denotes multiplication by g (x) and K N denotes the operator with kernel K N (x,y). 
This can be obtained by expanding out the product in the integrand, using the formula 

N\ r r 

^ N _ n y J---J Pn(xi, ...,x n , x n+1 , . . . , x N ) dx n+1 ■ ■ ■ dx N = det K N (x h xj) |™ j=1 



for the n-point correlation function, and then recognizing the resulting sum of multiple 
integrals as the expansion of the Fredholm determinant. Or it can be obtained by a simpler 
algebraic device 0. In our case 1 + g(x) = exp{f(2 1 ^ 2 N 1 ^ e (x — \f2N)/a)}. If we make the 
changes of variable 
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we find that the characteristic function equals the determinant of / plus the operator with 
kernel 



Now one has the scaling limit 

1 ^ / x /tttt y 



a; -y 

precisely the Airy kernel. Thus we see that the large N limit of the characteristic function 
equals 0(s) = det (J + A a (h)) where h(x) = e ls ^ x ^ — 1. Our asymptotic formula yields 

{isopl^ r°° s 2 r°° 1 

- J Q y/xf{-x)dx-— j x(g(x)) dx + o(l)\, 

where as before g(x) = f(—x 2 ). 

Notice that the limiting characteristic function is quadratic in s and hence the distribution 
is asymptotically normal. Of course this is not surprising since this occurs for other matrix 
ensembles and other scaling limits. Notice, though, that in this case the mean and variance 
of the limiting distribution only depend on the negative values of the argument of the original 
/. This is a reflection of the fact that the Airy function goes rapidly to zero for positive values 
and oscillates and tends to zero slowly for negative values of the argument. A question left 
to the future is how the asymptotics of functions A(x) and B(x) in a kernel of the form 

A(x)B(y) - A(y)B(x) 
x-y 

affect the asymptotics of the corresponding distribution functions. 
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